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ABSTRACT 

High resolution optical spectra of 57 Galactic B-type supergiant stars have been analyzed 
to determine their rotational and macroturbulent velocities. In addition, their atmospheric pa- 
rameters (effective temperature, surface gravity and microturbulent velocity) and surface ni- 
trogen abundances have been estimated using a non-LTE grid of model atmospheres. Com- 
parisons of the projected rotational velocities have been made with the predictions of stellar 
evolutionary models and in general good agreement was found. However for a small num- 
ber of targets, their observed rotational velocities were significantly larger than predicted, 
although their nitrogen abundances were consistent with the rest of the sample. We conclude 
that binarity may have played a role in generating their large rotational velocities. No corre- 
lation was found between nitrogen abundances and the current projected rotational velocities. 
However a correlation was found with the inferred projected rotational velocities of the main 
sequence precursors of our supergiant sample. This correlation is again in agreement with the 
predictions of single star evolutionary models that incorporate rotational mixing. The origin 
of the macroturbulent and microturbulent velocity fields is discussed and our results support 
previous theoretical studies that link the former to sub-photospheric convection and the latter 
to non-radial gravity mode oscillations. In addition, we have attempted to identify differential 
rotation in our most rapidly rotating targets. 
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early-type stars (see, for example, Maeder & Maynet 2001) and in 
long-duration gamma-ray bursts (Woosley & Heger uOOq) . In turn 
this h as stimulated several recent studies (Abt, Levato & Grosso 
2001 Huang & Gies |2004 Hun ter et al. l2008al: Martay an et al. 
20061 : Strom, Wolff & Dror 12001 Wolff et al l 19821 1200% of stel- 
lar rotation in both field and cluster stars and in different metallicity 
environments. However recently both slowly rotating stars with rel- 
atively high nitrogen abundance and rapidly rota ting stars wi th rel- 
atively low nitrogen abundances (Hunter et al. 2008b, 2009) have 
been identified in the VLT-FLAMES survey of massive stars (Evans 
et al. |2005). Hence other mechanisms as well as rotation may be 
important for mixing nucleosynthetically processed material to the 
surface, as discussed by Brott et al. d2008f) . 



In this paper, we present high resolution and high signal-to- 
noise spectra of a sample of 57 Galactic B-type supergiant stars, 
and use these data to further investigate the relationship between 
rotation and mixing, in addition to testing the predictions of theo- 
retical models. We also characterize the microturbulent and macro- 
turbulent velocity fields present in our sample and c ompare these 
with the theoretical studies of Cantiello et al. J2009I) and Aerts et 



1 INTRODUCTION 

The spectra of early-type stars are affected by a variety of non- 
thermal velocity fields. Those present in the photosphere are often 
desig nated as microturbulence (see, for example, Gies & Lambert 
1 1993 : Daflon, Cunha & Butler l2004 Simon-Diaz et al.l2006h and 



macr oturbul ence (see, for example, Howarth et al. 1 1997 ; Ryans 



et al. l2002h . Microturbulence is associated with distance scales 
smaller than the mean free path of a photon, while macroturbu- 
lence is associated with longer distance scales. Recently, the former 
has been linke d with the presence of sub-surface convection fields 
(Cantiello et al. 2009) driven by the opacity of iron group elements, 
whilst the latter may be a manifestation of the large number of non- 
radial gravity-mode stellar oscillations (Aerts et al. 1 20091) present in 
the stellar photosphere. 

Stellar rotation also significantly affects early-type stellar 
spectra and particularly for main sequence objects can dominate 
the broadening o f the metal absorption lines (see, for example, 
Grav ll97fj] 2008). Rotation may play a key role in the evolution of 
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al. (§009) respectively. We have also attempted to make the first 
identification of differential rotation in evolved massive stars. 



2 OBSERVATIONS AND DATA REDUCTION 

The observing list was initially developed to sample as widely 
as possible B-type supergiants of different spectral subtype. The 
availibility of targets (in terms of spatial position and brightness) 
led to the early-B spectral types being better sampled than the 
later spectral types as can be seen from TableQ] However even 
for the latter, most spectral subclasses contained at least one tar- 
get. A number of targets were subsequently excluded from the 
analysis. In particular, four B9 supergiants (HD 94367, HD 96919, 
HD 1 1 1904, HD 158799) were found to have effective temperatures 
that were smaller than the lower limit of our model atmosphere 
grid. Additionally eight targets (HD 74804, HD 99103, HD 75821, 
HD76728, HD 101330, HD 152046, HD 68161, HD 53244) were 
found to have gravities that were significantly higher than supe- 
giants of similar spectral types. Indeed their logarithmic gravity es- 
timates (3.0-3.5 dex) were similar to those found for B-type giants 
(see, for example, Vrancken et al. 1200 0) and they were therefore 
excluded. 

Previous investigations of B-type supergiants (Ryans et al. 

120021 ; Dufton et al. l2006l : Simon-Diaz & Herrero l2007h have shown 
that many have relatively small projected rotational velocities. In 
order to try to increase our sample of rapidly rotating supergiants 
(in order to, for example, investigate whether they have a differ- 
ent evolutionary history), we have used the projected rotational ve- 
locities estimated by Howarth et al. ( 1997) from IUE spectra, as 
a guide in selecting additional targets. These were determined by a 
cross-correlation technique, which implicity assumed that the spec- 
tral line broadening was dominated by rotation. Howarth et al. dis- 
cussed the validity of this assumption and inferred from the lack of 
supergiants with very small projected rotational velocities that an- 
other important broadening mechanism was probably present. In- 
deed for the majority of their B-type supergiants this mechanism 
probably dominates as illustrated by Fig. Q] which plots their esti- 
mated projected rotational velocity against spectral type. Most of 
the sample lie on a locus with an upper bound stretching from ap- 
proximately 100 km s _1 at a B0 spectral type to 75 km s _1 at 
B3. For these targets the projected rotational velocities are prob- 
ably small and the broadening is dominated by macroturbulence. 
This is supported by fig. 2 of Dufton et al. (2006), where a similar 
dependence of the macroturbulent broadening with spectral type 
was found. However several supergiants have estimates that appear 
to be too large for their spectral type and in these cases the total line 
broadening is probably being significantly augmented by rotation. 
Four of these targets had already been included in our observing 
lists but six additional targets had positions and apparant magni- 
tudes that made them accessible and these were observed and have 
been identified with an asterisk in the TableQ] 

High resolution (R~48,000) spectra were obtained with the 
Fibre fed Extended Range Optical Spectrograph (FEROS; Kaufer 
et al. 1 19991) on the 2.2-m MPI/ESO telescope at La Silla during a 
three night observing run in April 2005. Some preliminary observa- 
tions had been obtained in July and December 2004 to test the fea- 
sibility of our analysis procedures (see Sect. [3} but due to their rela- 
tively low signal-to-noise-ratios (SNRs) these have not been used in 
the current analysis. All spectra covered the wavelength range from 
3600 A to 9200 A and were reduced using the reduction pipeline 
that runs under the MIDAS environment (Kaufer et al. 1999). 
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Figure 1. Estimated projected rotational velocities as a function of spectral 
type for the early B-type supergiants analysed by Howarth et al. 419971) . The 
dashed line represents a subjective estimate of the upper limit for the part 
of the sample that are rotating relatively slowly. 



The signal-to-noise ratio (SNR) for each spectrum was esti- 
mated by taking a normalised region of the continuum and finding 
the standard deviation, the SNR estimate being its reciprocal. This 
procedure was normally carried out in the region of 4500-4600 A 
as this was compatible with the spectral lines used to estimate the 
projected rotational velocities (see Table |2). However for spectral 
types later than B5, this region contained a significant number of 
strong lines and hence the region from 4700-4800 A was adopted. 
As the SNR is expected to vary with wavelength, the estimates for 
several stars were measured in both regions with those from the 
4700-4800 A region being found to be approximately 5 per cent 
lower. Hence the estimates for all stars of spectral type B6 or later 
have been scaled accordingly. Given that weak undetected absorp- 
tion will lead to underestimates of the SNRs, these estimates should 
only be taken as an approximate comparative guide to the quality of 
the spectra. All the targets are listed in the TableQ] together with ex- 
posure times and estimates of the SNRs. Spectral types are taken in 
the first instance from the Bright Star Catalog (Hoffleir. ll99ll) . then 
from the IUE Atlas of B-Type Stellar Spectra (Walborn, Parker & 
Nichols. Il995h , the Uni versit y of Michigan Catalogue of Spectral 
Types (Houl & Cowley. [19751) and Reed ( 120031) . 

For all stars in the sample, the equivalent widths (Wj) of se- 
lected lines were estimated. The N II line at 3995 A was measured 
to estimate nitrogen abundances, whilst the silicon spectrum (Si II 
4128 A, 4131 A; Si III 4552 A, 4567 A, 4574 A; and Si IV 4116 A) 
was used to constrain effective temperatures as described in Sec- 
tion|4] The continuum adjacent to each line was normalised using a 
second-order polynomial fit. For the slowly rotating stars (see Sect. 
13. Il l where macroturbulent broadening will generally dominate, a 
Gaussian was then fitted to the absorption line using a least squares 
technique. For those stars with a relatively large projected rotational 
velocity (~ 100 km s _1 or greater), the lineshape is dominated by 
rotation. In these cases, a profile shape appropriate to rotationally 
dominated broadening (with the value of v e sin i taken from Table 
[3]- see Sect |3. It was used, and was generally found to give a bet- 
ter fit. Although the choice of profile is a possible source of error, 
we note that the equivalent width estimates of the stronger lines 
normally agreed to within 5 to 10 mA with those found by direct 
numerical integration over the profile. The estimates are presented 
in Tableg] 



Rotational velocities of Galactic B-type supergiants 3 



Table 1. Observational details. 



Identifier 


Sp. Type 


Exp. (s) 


SNR 


HD 77581* 


BO la 


716 


430 


HD 122879 


BO la 


464 


510 


HD 149038 


BO la 


119 


470 


HD 167264 


BO lab 


178 


500 


HD 164402 


B0 lb 


258 


550 


HD 168021 


B0 lb 


684 


400 


HD 115842 


B0.5 la 


321 


480 


HD 150898 


B0.5 la 


214 


500 


HD 152234 


B0.5 la 


190 


450 


HD 152667* 


B0.5 la 


362 


350 


HD 94493* 


B0.5 lab 


1010 


520 


HD 64760 


B0.5 lb 


62 


480 


HD 103779 


B0.5II 


962 


500 


HD 155985 


B0.7Ib 


486 


480 


HD 109867 


Bl la 


397 


460 


HD 152235 


Bl la 


424 


420 


HD 142758 


Bl la 


830 


410 


HD 148688 


Bl Iae 


170 


380 


HD 154090 


Bl la 


111 


470 


HD 152236 


Bl Iape 


98 


440 


HD 119646 


Bl Ib-II 


543 


410 


HD 125545 


Bl lab-lb 


1145 


400 


HD 150168 


Bl lab-lb 


229 


390 


HD 99857* 


Bl lb 


1220 


400 


HD 157246 


Bl lb 


27 


560 


HD 93840* 


BN1 lb 


1640 


410 


HD 106343 


B1.5 la 


386 


400 


HD 148379 


B1.5 Iape 


173 


400 


HD 96248 


BC1.5 lab 


528 


420 


HD 108002 


B2 la-lab 


750 


400 


HD 52089 


B2II 


5 


430 


HD 99953 


B2 lab-lb 


495 


440 


HD 93827* 


B2 lb -II 


6770 


270 


HD 111990 


B2Ib 


647 


400 


HD 117024 


B2Ib 


861 


400 


HD 165024 


B2Ib 


37 


510 


HD 141318 


B2II 


246 


450 


HD 92964 


B2.5 Iae 


178 


440 


HD 116084 


B2.5Ib 


270 


500 


HD75149 


B3Ia 


192 


440 


HD53138 


B3Iab 


20 


500 


HD51309 


B3II 


70 


520 


HD 157038 


B4 la 


460 


450 


HD 1591 10 


B4Ib 


1045 


550 


HD 79186 


B5 la 


126 


440 


HD 111973 


B5Ia 


288 


480 


HD 58350 


B5Ia 


12 


440 


HD 86440 


B5Ib 


33 


520 


HD 164353 


B5Ib 


49 


520 


HD 83183 


B5 II 


54 


400 


HD 74371 


B6Iae 


155 


440 


HD 80558 


B6Iae 


280 


430 


HD 105071 


B6 lab-lb 


428 


510 


HD 125288 


B6 lb 


68 


500 


HD91619 


B7Iae 


359 


460 


HD 111558 


B7Ib 


1007 


450 


HD 166937 


B8 lap 


44 


460 



4+ 



3 
It 



* Subsequently added to observing list as Howarth et al. (1997) 

found anomalously large line widths for spectral type 

t Listed as binary in SIMBAD online database. 

(1) Hoffleit & Warren (1991) . (2) Walborn, Parker & Nichols (1995) 

(3) Houk & Cowley (1975). (4) Reed (2003) 



3 DATA ANALYSIS 

3.1 Projected rotational velocities 

The use of Fourier transforms to estimate projected rotational ve- 
locit ies {v E sin i) from stellar spectra was first proposed b y Car- 
roll dl933h . and more recently has been dis cussed by Gray Gaze, 
I2008I) . Reiners & Schmitt d2002h and Royer d2005h . The applicabil- 
ity of this methodology speci fically to OB-type stars has been con- 
sidered by Simon-Diaz et al. ( 120061) . In essence the Fourier Trans- 
form technique is predicated on the convolution theorem, whereby 
when an observed spectrum is transformed into the Fourier domain, 
the convolution of the intrinsic spectrum with the rotational, macro- 
turbulent and instrumental profiles becomes a multiplication of the 
corresponding Fourier transforms. Furthermore, as only the rota- 
tional profile is expected to have zeroes in its Fourier Transform at 
low frequencies, these will appear in the total transform, thereby 
allowing the projected rotational velocity to be estimated. 

The position of the first zero for a solidly rotating body is in- 
versely prop ortional to th e projected rotational velocity (see, for ex- 
ample, Gray 1976, 200il); this relationship will form the basis for 
our estimates of the projected rotational velocities. A more sophisti- 
cated analysis, such as that discussed by Reiners (2003) shows that 
differential rotation, limb darkening, the angle of inclination and 
gravity darkening can all alter the line profile, and hence its Fourier 
transform. However, these effects should not significantly affect the 
position of the first zero, particularly for stars with an equatorial ve- 
locity, v < 200 km s _1 which is appropriate to most of our sam- 
ple. Furthermore, the main diagnostic of these secondary effects - 
namely the ratio of the first and second zeroes in the Fourier trans- 
form, (qi /q2 ) - was only measurable for stars which had significant 
rotational broadening. Such objects will have large equatorial ve- 
locities and are discussed further in Sect. |5.6l 

As discussed in Sect.[2] there is additional broadening present 
in the spectr um of early-type supergia nts (H owarth et al. 1 19971 ; 
Ryans et. al. 120021 ; Simon-Diaz et. al. 120061) . which is normally 
designated as macro turbulence. This will also affect the line pro- 
file in the Fourier domain, acting to decrease the amplitude of the 
sidelobes. While the precise nature of this broadening is not fully 
understood, if it can be approximated by, for example, a Gaus- 
sian it should not affect the postio n of th e first zero in the Fourier 
Transform. Recently Aerts et al. j2009h simulated the effects of 
non-radial gravity mode oscillations on the observed line profile. 
In particular, they find that on occasions this can lead to an under- 
estimation of the projected rotational velocity and we discuss this 
further in Sect. [5] 

Thirteen metal spectral lines were considered for each star, as 
listed in Table [2] Four of the lines used are in close doublets or 
triplets with separations of 0.25 A or less (corresponding to a ve- 
locity of approximately 15 km s _1 ), which is less than the typical 
broadening from macroturbulence and rotation. Indeed, in all these 
cases, the features appeared as single lines. Hydrogen and diffuse 
helium lines, although strong and well observed, were not normally 
considered, as although the Fourier technique should be capable of 
separating rotational broadening from that due to linear Stark ef- 
fect, the approximation of line broadening as a convolution of unre- 
lated effects may become unreliable (Heinzel. ll978l) . An exception 
was made for two targets (HD 64760 and HD 152667) with large 
projected rotational velocities, where the metal lines were poorly 
observed and/or blended; for these targets, the diffuse neutral he- 
lium lines at 3819A and 4026A were also considered. 

Spectra were analysed using PROCSPEC, an unpublished pack- 
age of IDL routines for manipulating spectra. The lines listed in Ta- 
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ble [2] were normalised using a polynomial that had been fitted to 
the adjacent continuum. Whilst the amount of continuum consid- 
ered should be minimised to reduce the noise in the Fourier do- 
main, tests showed that the projected rotational velocity estimates 
were not significantly affected by this factor. 

Each line was then assigned a ranking between 1 and 5 to re- 
flect its reliability (relative to the other lines for that star). Factors 
considered were the quality of the spectrum, degree of symmetry in 
the line, the presence of cosmic rays, confidence in the normalisa- 
tion procedure, and the intrinsic strength of the line. Lines that were 
in the top two quality ranking bins for each star were Fourier trans- 
formed using the procedure discussed by Simon-Diaz & Herrero 
(2007). Those that gave a transform with a clearly identifiable first 
zero were then used to estimate the projected rotational velocity. 
An illustrative example of this procedure can be seen in Fig. [2] For 
each star, the number, mean and standard deviations of these esti- 
mates of the projected rotational velocity are summarized in Table 
® 



Table 2. Lines used to estimate the projected rotational velocities and 
macroturbulence for our targets. 



Species Wavelength (A) 



Nil 


3995.00 


Sill 


4128.07 


Sin 


4130.89 


CII 


4267.00 / 4267 


Mgll 


4481.13/4481 


Sim 


4552.62 


Sim 


4567.82 


Sim 


4574.76 


on 


4590.97 


on 


4595.96/4596 


Nil 


4630.54 


on 


4661.63 


He I 


4713.15/4713 



3.2 Macroturbulence 

Macroturbulence was originally introduced as an additional veloc- 
ity field to improve the agreement between observed metal line pro- 
files of supergiants and those calc ulated using model atmos pheres 
(see, for ex ample , Howarth et al. 1 19971 ; Simon-Diaz et al. 2006; 
Ryans et al. l2002ri . Macroturbulence is postulated to be a velocity 
field characterised by a length scale longer than the mean free path 
of a photon (as opposed to microturbulence). Recent work by Aerts 
et al. (2009) has suggested a physical explanation for macrotur- 
bulence, in term of the pulsation of stars due to non-radial gravity- 
mode oscillations. These authors considered pulsational models for 
a Bl-type supergiant, with comparable atmospheric parameters to 
those found in this paper and found that the additional line broad- 
ening simulated a macroturbulence consistent with those found ob- 
servationally. Additionally the pulsations could lead to errors in the 
values estimated for projected rotational velocities and we return to 
this in Sect.|5H 

We have estimated macroturbulences for all the stars in our 
sample by assuming that the magnitude of the velocity field fol- 
lows a Gaussian distribution. Clearly given the work of Aerts et al., 
this assumption may not be valid. However the estimates will still 
be useful in characterising how the excess broadening varies with 
spectral type and in many cases where macroturbulence dominates 
the line broadening, the observed profiles appeared to be well ap- 
proximated by a Gaussian function. 

Our procedure was to consider the lines tabulated in Table [2] 
(excluding the neutral helium line) and to optimise the agreement 
between theoretical and observed profiles. Intrinsic spectra were 
extracted from our TLUSTY grid using the grid point closest to our 
estimated atmospheric parameters and an elemental abundance that 
gave an equivalent width similar to that observed - typically differ- 
ences were less than ten per cent. The theoretical profile was then 
scaled so that the two equivalent widths were the same. The theo- 
retical profile was convolved with a rotational broadening function 
using the values of the projected rotational velocity listed in Table 
[3]and then with Gaussian profiles to represent the macroturbulence. 
For each line an estimate of the macroturbulence was found by min- 
imising the sum of the square of the residuals with the mean, stan- 
dard deviation and number of the estimates being listed in Table 
[3] Note that we have defined the macroturbulence as the -| (half-) 
width of the Gaussian profile. These values are discussed further in 
Sect.El 
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Figure 2. An example of the Fourier transform methodology for the CII 
4267A line in HD1 17024. The topmost panel shows the observed linepro- 
file, normalised to the continuum. The centre panel shows the Fourier trans- 
form of this line. The bottom panel shows the same FT, but with two theoret- 
ical Fourier transforms overplotted, one with no macroturbulence (red line), 
and one with a Gaussian macroturbulent velocity of 50km s _1 (blue line). 
The value of v e sin i which gave the best fit to this line was 53 km s — 1 . 
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Table 3. Estimates of the projected rotational velocity, atmospheric parameters and nitrogen abundances. The ordering of targets follows that in Tabl^T] 



Identifier 


v e sin i 




n r 


Vf 






T e ff 






r— i 


M 




km s _1 


km s — 1 




km s -1 


km s _1 




K 


cm s -2 


km s — 1 


M Q 


un 77S8 1 ++ 


JO 


1 


4 








26500 


2 90 


17 


8 40 


40** 


RD 1 9987Q 


J J 


a 

j 


u 


oU 


7 


5 


27200 


9 QO 


14 


8 37 

O.J / 


40** 


un 1 4Q018 


u 1 


a 


J 


6Q 


8 


j 


28500 


1 10 
J.1U 


12 


7 Q1 


40** 


RD 1 67964 


8^ 


7 




66 


7 


5 


27500 


3 05 


16 


7 96 


30** 


RD 1 64409 


46 


5 


7 


56 


7 


7 


28100 


3.25 


13 


7 55 


30** 


RD 1 68091 ++ 

lily 1UOUZ.1 


48 


9 


7 


58 


7 


7 


26900 


3.10 


15 


7 81 


25** 


RD 1 1 S849 

111./ 1 IJOti 


39 




5 


66 


5 


5 


24800 


2.75 


14 


8.44 


40 


RD 1 S08Q8 

lily 1JU070 


93 


4 


4 


56 


9 


5 


26400 


3.10 


20 


8.00 


25** 


RD 1 59914++ 

111./ IjAiJt 


59 


3 


4 


57 


3 


8 


25700 


2.90 


14 


7 69 


38** 


RD 1 S9667 
ni/ ijiuu / 


1 39 




<^ 








76000t 


3 05 


18 


7 58 


25** 


HD 94493 


97 




7 


60 


5 


g 


23500 


2 95 


14 


8 02 


24** 


rd 64760 


-jj 


7 










9 6000 t 


1 9S 


22 


7 89 


91 

L J 


RD 1 (Y377Q 
j 1 1./ i ij j i i y 


40 


1 




60 


4 


5 


25000 


3.00 


12 


8 07 


25** 


RD 1 SSQ8S 
ni./ i JJ70J 


35 


4 


9 


46 


3 




23200 


2 95 


12 


8 19 


22 


RD 1 09867 


43 




9 


56 


4 


9 


22300 


2.75 


14 


8.32 


23 


un 1 59915 


46 




7 


58 


7 


8 


22000 


2.65 


14 


7 59 


30 


un I497S8 
ni.' itz, / jo 


46 


g 




18 

JO 


7 




17300 


9 9S 


21 


8 46 


is 

j j 


RD 148688 


48 


5 


7 


44 


5 


9 


20700 


2 45 


16 


7 94 


37 


un 1 S40Q0 


47 


s 

j 


4 


54 


4 


□ 


22200 


2 70 




8 10 


"^8 


RD 1 S9916 




a 

j 




46 







21500* 










un 1 i q^4A 


jj 






96 


4 


Q 



17600 


2 50 


2 j 


7 Ql 

/ .7J 


20 


RD 1 9SS4S 


65 




7 


53 


4 


9 


90800 


2 75 


17 


7 70 


22 


un 1 snifis 

\ \ \J 1JU1UO 


129 


7 




55 


12 


5 


24800 


3 15 


1 3 


8.26 


20** 


HD 99857 

i ny yyoji 


189 


10 


6 


59 


1 j 


7 


21500 t 


2.95 


16 


7.20 


18 


RD 1 S7946 


269 


1 1 


4 








21500^ 


2.90 


15 


8.1 1 


19 


HD 93840 

Illy y^jyj^xj 


58 


3 


6 


68 


4 


6 


20900 


2.75 


15 


8.37 


23 


un 1 06141 

lily i uujtj 


43 


4 




52 


4 


7 


20 1 00 


2.50 


17 


8.20 


30 


HD 148379 

Illy IIOJ / 7 


44 






33 


4 


8 


17000 


2.00 


19 


7.93 


43 


HD 96748 
i \ \y 7UZito 


i7 


5 


7 


5 3 


4 


g 


19500 


2 40 


15 


7 48 


29 


hd 1 osoo? 


37 


4 


8, 


46 


4 




20200 


2.60 


16 


8.19 


25 


un 5908Q 
nu Jiuo? 


22 


3 




17 




g 


20100 


3 05 


17 


7 74 


14 


1 1 1 .V yyyjj 


49 




5 


37 


4 


5 


16800 


2 15 


22 


8 43 


32 


RDQ1897 
1 1 iy 7joi/ 


227 


14 


4 








18500"t 


2 70 


20 


7 86 


17 


RD 1 1 1 QQO 
1 1 iy iii 77U 


34 


6 


7 


30 


g 


9 


16500 
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t large u e sin i precludes use of silicon ionization balance and u c sin i estimates supplemented with those from He lines 

* gravity is too low to be analysed with TLUSTY grid 

* * In ambiguous region of HR diagram that can lead to ambiguous mass estimates 
+ Si III lines too weak to reliably estimate the microturbulence 

++ Classified as binary in SIMBAD 
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Table 4. Equivalent widths (in mA) for selected metal lines. The ordering of targets follows Table[T] The letter B indicates a blend of lines. 



Identifier 


N II 3995 


Si IV 41 16 Si II 4129 


Si II 4131 


Si III 4553 


Si HI 4568 


Si III 4575 


HD 77581 


155 


475 




435 


365 


195 


HD 122879 


105 


455 




305 


230 


115 


HD 149038 


60 


460 




230 


175 


85 


HD 167264 


75 


460 




270 


200 


95 


HD 164402 


50 


345 




270 


210 


110 


HD 168021 


85 


365 




300 


245 


125 


HD 115842 


165 


395 




405 


345 


195 


HD 150898 


155 


370 




325 


250 


125 


HD 152234 


80 


380 




350 


295 


160 


HD 152667 


85 


B 




495 


400 


230 


HD 94493 


205 


190 




375 


310 


180 


HD 64760 


170 


B 




410 


330 


155 


HD 103779 


165 


255 




345 


285 


165 


HD 155985 


235 


165 




375 


315 


195 


HD 109867 


280 


180 




425 


355 


215 


HD 152235 


125 


230 




495 


425 


270 


HD 142758 


385 


85 




555 


470 


305 


HD 148688 


210 


200 




550 


480 


305 


HD 154090 


285 


190 




470 


400 


245 


HD 152236 


470 


95 




500 


435 


275 


HD 119646 


235 


40 




380 


315 


185 


HD 125545 


200 


135 




505 


425 


270 


HD 150168 


230 


230 




385 


330 


195 


HD 99857 


95 


B 




405 


330 


200 


HD 157246 


290 


B 




460 


375 


255 


HD 93840 


350 


105 




420 


345 


215 


HD 106343 


330 


115 




465 


385 


235 


HD 148379 


240 


40 




430 


365 


230 


HD 96248 


140 


95 




465 


395 


240 


HD 108002 


315 


115 




450 


380 


235 


HD 52089 


200 


55 




350 


295 


190 


HD 99953 


375 


60 




460 


380 


235 


HD 93827 


225 


B 




390 


325 


210 


HD 111990 


215 


20 




275 


215 


130 


HD 117024 


130 


20 




330 


270 


160 


HD 165024 


300 


55 




405 


340 


215 


HD 141318 


220 


30 




310 


250 


160 


HD 92964 


345 


25 




365 


300 


180 


HD 116084 


240 


30 




370 


305 


180 


HD75149 


285 


140 


170 


255 


200 


115 


HD 53138 


200 


155 


180 


230 


175 


105 


HD 51309 


120 


150 


160 


165 


125 


70 


HD 157038 


340 


170 


190 


300 


250 


155 


HD 159110 


95 


65 


65 


195 


150 


90 


HD 79186 


195 


155 


170 


210 


160 


95 


HD 111973 


190 


150 


170 


230 


180 


105 


HD 58350 


150 


190 


195 


155 


115 


65 


HD 86440 


60 


165 


175 


85 


55 


25 


HD 164353 


100 


175 


185 


130 


90 


50 


HD 83183 


55 


175 


185 


75 


45 


20 


HD 74371 


80 


205 


215 


115 


80 


45 


HD 80558 


105 


200 


205 


95 


65 


35 


HD 105071 


55 


250 


255 


65 


35 


15 


HD 125288 


50 


170 


175 


60 


30 


15 


HD91619 


135 


200 


210 


110 


80 


45 


HD 111558 


55 


305 


315 


50 


30 


15 


HD 166937 


50 


250 


260 


60 


40 


20 
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4 ATMOSPHERIC PARAMETERS 

Atmospheric parameters and nitrogen abundances have been esti- 
mated from non-LTE model atmosphere grids generated usi ng the 
codes TLUST Y and S YNSPEC (Hubeny 1988; Hubeny & Lanz l 19951 : 
Hubeny et al. 1 19981 : Lanz and Hubeny 2007). Details of the meth- 
ods adopted can be found in Hunter et al. l2007h while a more de- 
tailed discussi on of t he grids can be found in Ryans et al. (l2003h and 
Dufton et al. (2005f|. Hence only a brief summary will be given 
here with further details available as discussed above. 

Four model atmosphere grids have been calculated for metal- 
licities correponding to a Galactic metallicity of [Fe/H]=7.5 dex, 
and metallicities of 7.2, 6.8 and 6.4 dex to represent the LMC, SMC 
and lower metallicity material respectively. For each of these grids 
non-LTE models have been calculated for effective temperatures 
ranging from 12 000K to 35 000K, in steps of no more than 2 500K, 
surface gravities ranging from 4.5 dex down to the Eddington limit, 
in steps of no more than 0.25 dex and for microturbulences of 0, 
5, 10, 20, and 30kms _1 . Assuming that the light elements (C, 
N, O, Mg and Si) have a negligible effect on line-blanketing and 
the structure of the stellar atmosphere, models were then generated 
with the light element abundance varied by +0.8, +0.4, -0.4 and - 
0.8 dex about their normal abundance at each point on the TLUSTY 
grid. Theoretical spectra and equivalent widths were then calcu- 
lated based on these models. Photospheric abundance estimates for 
approximately 200 absorption lines at any set of atmospheric pa- 
rameters covered by the grid can then be calculated by interpolation 
between the models via simple IDL routines. The reliability of the 
interpolation technique has been verified by Ryans et al. ( 120031) . 

As the estimation of the atmospheric parameters is inter- 
related, it was necessary to use an iterative approach, as follows: 

- initial estimates of the effective temperature (T c a) and gravity 
(log g) were obtained from the spectral type and the calibration of 
Crowther, Lennon & Walborn J2006l) . 

- the microturbulence (£ t ) was estimated from the Sim triplet at 
4552-4574A. 

- the effective temperature was re-estimated from the silicon 
ionization equilibrium. 

- surface gravity was re-estimated from the Balmer line profiles. 

- the last three steps were iterated until the estimates converged. 

Further details of this approach are given below. 



4.1 Effective temperature 

For the hotter targets with spectral types earlier than B3, the Si III 
to Si IV ionization balance was used, whilst for the later spectral 
types the Si II to Si III ionization balance was considered. The very 
high quality of the observational data implies that the random er- 
rors should be small and normally less than 1 000 K. Systematic 
errors due to, for example, errors in the adopted atomic data or 
in the physical assumptions were more difficult to quantify. How- 
ever, Dufton et al. J2005h compared atmospheric parameters found 
for SMC supergiants using the approach adopted here with those 
found for two stars by Trundle et al. (2004). The latter utilized the 



unified code FASTWIND ( Santo laya-Rey, Puis & Herrero, A. 19971 : 
Herrero, Puis & Naiarro I2OO2I : Repolust, Puis & Herrero l2004h . 
The two approaches yielded similar estimates of the effective tem- 
perature (agreeing to 500K for one star and 2000K for the other), 
hence an error estimate of ± 1 000 K would appear to be appropri- 
ate. The values of the other atmospheric parameters from the two 
codes were also in good agreement; log g to within 0.1 dex, micro- 
turbulence to within 3 km s _1 and nitrogen abundance to within 0.2 
dex. 

Five targets (HD 152667, HD 64760, HD 99857, HD 157246, 
HD 93827) had large projected rotational velocities that precluded 
the observation of two ionization stages of silicon. In these cases, 
the effective temperature was taken to be that implied by the spec- 
tral type and then the other atmospheric parameters were estimated 
as discussed below. 



4.2 Surface Gravity 

The logarithmic surface gravity (log g) of each star was estimated 
by fitting the observed hydrogen Balmer lines with theoretical pro- 
files. In order to minimise the effects of the stellar wind, relative 
high order lines in the series (H7 and H<5) were normally consid- 
ered. Automated procedures have been developed to fit model spec- 
tra in our TLUSTY model atmosphere grid to the observed spec- 
tra, with contour maps displaying the region of best fit. With an 
effective temperature estimate being available (from the methods 
described above), it was a straightforward matter to estimate the 

gravity. 

Hunter et al. J2007l) have discussed the errors associated with 
this procedure and estimate a typical random error of ±0.1 dex. 
Given that our observational data is generally of a higher quality 
than their Magellanic Cloud spectroscopy, this would appear to be 
a conservative error estimate to adopt here. Systematic errors due 
to, for example, uncertainties in the line broadening or assumption 
in the model atmosphere calculations could be significant, although 
the gravity estimates from the TL USTY grid and FASTWIND calcu- 
lations discussed in Dufton et al. (2005J) agreed to within 0. 1 dex. 
Hence an error estimate of ±0.2 dex would therefore appear to be 
a conservative estimate. 



4.3 Microturbulence 

The microturbulence is normally derived by removing any system- 
atic dependence of the abundance estimates on line strength found 
from lines of a specific ion. For B-type stars, the O il ion is of- 
ten co nsidered (for ex ample , Simon- Diaz et al. 2006; Hunter et al. 
120051 : Gies & Lambert fl992 : and Daflon et al. l2004h as its rich spec- 
trum should improve its reliablity. However the analysis is com- 
plicated by the lines arising from different multiplets, making any 
estimate susceptible to errors in the adopted atomic data or in the 
magnitude of non-LTE effects. In order to remove these uncertain- 
ties, a single Oil multiplet can be used. However the number of 
lines is then drastically reduced, whilst the O II spectrum is weak 
or unobservable for spectral types later than B3. 

In order to maintain consistency throughout the analysis we 
have instead estimated the microturbulence from the Si III triplet 
of lines at 4552^1-574 A, which is observed in all our spectra. As 
these lines are from the same multiplet, errors arising from the os- 
cillator strengths and non-LTE effects should be negligible and this 



See also http://star.pst.qub.ac.uk 



metho d has been used pr eviously by, for example, Dufton et al. 
( 2005) and Vrancken et al. dl997li2000h . For two stars (HD 10507 1, 
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HD 125288), the Si III lines were too weak to usefully constrain the 
microturbulence. For these cases a value of 20 km s~ 1 was adopted, 
which is compatible with the estimates for other supergiants of sim- 
ilar spectral type. 

4.4 Nitrogen abundances 

In early B-type stars, there is a relatively rich N II spectra with sin- 
glet transitions at 3995 and 4447 A and triplet multiplets at 4601- 
4643, 4780-4819 and 4994-5007 A. Here we have estimated nitro- 
gen abundances solely from the feature at 3995 A. This choice was 
based on the following considerations: 

(i) this feature is amongst the strongest N II lines in the optical 
spectrum and is unblended. As such it could be measured in effec- 
tively all our targets (see Tabled- 

(ii) the other lines and multiplets were either not observable over 
the full range of spectral and/or suffered from blending. 

(iii) the use of only one feature will inevitably increase the ran- 
dom error in the estimate and may lead to systematic errors. In the 
case of the former, the abundance estimates listed in the Table [3] 
have a range of over 1 dex, which is larger than our estimated ran- 
dom error discussed below. Additionally the nitrogen abundances 
will depend on the amount of nucleosynthetically processed ma- 
terial mixed from the core to the photosphere. Hence the relative 
rather than the absolute abundances will be important for investi- 
gating this phenomena. 

The nitrogen abundance estimates form the line at 3995 A 
were obtained using our grid of TLUSTY models and the methodol- 
ogy discussed, for example, in Hunter et al. fe007h . 

To investigate the magnitude of the errors in these abundances, 
we have also deduced nitrogen abundance from the other N II lines 
listed above for six stars covering the range of spectral types (and 
hence effective temperatures present in our sample). These are sum- 
marized in Table [5] where to ease comparison we have also listed 
the estimates obtained from the N II line at 3995A, which were 
taken directly from Table [3] As can be seen there is no system- 
atic differences between the two set of estimates with the differ- 
ences being compatible with the standard deviations found from 
the abundances estimates from the other features. Hence we do not 
believe that our choice of a single feature has introduced significant 
systematic errors. The standard deviations in Table [5] range from 
±0.1 1-0. 17 dex and will represent at least in part the random errors 
in the abundance determinations. Also there may be systematic dif- 
ferences in the estimates from different multiplets which will con- 
tribute to the standard deviations of the corresponding estimates. 
Hence we adopt a random error of ±0.15 for the estimates in Ta- 
ble [3] We believe that this may be conservative as the feature at 
3995 A is stronger and does not suffer from the blending affecting 
some other feature. As well as random errors there could addition- 
ally systematic errors due to for example errors in the atmospheric 
param eters. These have been extensively discussed by Hunter et al. 
(120071) and vary depending on the strength of the feature and on 
the adopted atmospheric parameters. Typically such errors are of 
the order of 0.1-0.2 dex (but on occasions can be larger). Here we 
will adopt a typical error in our nitrogen abundance estimates of the 
order of 0.2 dex. 

4.5 Masses 

Masses for the stars in our sample were estimated by comparing our 
effective temperatures and surface gravities estimates with those 



Table 5. Nitrogen abundances estimated form the feature at 3995 A (N3ggs ) 
and from other N II lines (N ) for six targets covering the range of effective 
temperatures of our sample. For the latter the standard deviation (cr) and 
number (n) of the individual estimates are also listed. 
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(7 
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0.17 


8 
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12 
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17300 


8.46 
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HD 79186 


15100 


8.39 


8.43 
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Effective Temperature (log K) 
Figure 3. Estimated effective temperatures and surface gravities for our 
sample, plotted against the Geneva group evolutionary tracks (Meynet & 
Mader, 2003). Each track is identified by an integer which is the initial 
mass in solar masses (Mq). The number in brackets refers to the mass of 
the model at an effective temperature T c s= 4.1, note that the mass of the 
models changes by less than 0.2 Mq between the loop back at the upper 
left and this point. 

predicted from the evolu tonary tracks from the Geneva group mod- 
els (Meynet & Maeder, 120031) . All the stars in our sample were 
found to lie in a region between approximately IOMq and 40Mq. 

As can be seen from Fig. [3] the evolutionary tracks in our 
pseudo-HR diagram loop back on themselves for values of log T c s 
greater than 4.3 dex and log g greater than 2.8 dex. This increases 
the uncertainty in the determination of mass for stars in this region 
of the diagram, as a degeneracy in log g vs T a ff is introduced. This 
uncertainty is mitigated by the fact that the stars in this region all 
have masses between 20Mq and 40Mq and the behaviour of evo- 
lutionary models in this mass range is quite similar. 

At higher temperatures and lower gravities, the tracks are 
closer together, leading to larger errors in the mass estimates, with 
the relative error being approximately fifteen per cent. Compound- 
ing this problem, stars with high rotational velocities are likely to 
exhibit grav ity darkening in accordance with Von Zeipel's Theorem 
(Von Zeipel [l924al 1 1 924bh . and if these stars are pole on, signifi- 
cant err ors can be made in estimating their effective temperature 
(Gillich, 2008). These uncertainties in effective temperature will in 
turn increase the uncertainty in mass. Based on these factors, we 
have adopted a conservative uncertainty in mass of ± 20 per cent, 
which corresponds to an uncertainty of 2Mq at the lower end of 
our mass range and 8M© at the upper end. 
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5 RESULTS AND DISCUSSION 

5.1 Comparison with previous rotational velocity estimates 

There are several studies of projected rotational velocities that have 
target s in co mmo n with our sample (see, for e xample, Abt et al. 
|2002| ; Pennv ll996l ; Simon-Diaz & Herrero [2007h . It is encouraging 
that Simon-Diaz & Herrero deduced from a similar methodology 
to that used here a value of 73 ± 4 km s _1 for HD 167264. This is 
in reasonable agreement with our estimate of 83 ± 3 km s (our 
error estimate for the mean is based on the standard deviation of 
the sample and assuming that the errors are normally distributed). 
Kaufer, Prinja & Stahl ( 120021) have found v c sin i to be 265 ± 5 
km s^for HD 64760, which is again in good agreement with our 
measured value of 255 km s~ . It is also encouraging that Kaufer et 
al. did not use the Fourier method to obtain this result, but rather ro- 
tationally broadened a model spectral line to match the Sim line at 
4552 A. Unfortunately we had no targets in common with the sam- 
ple of Markova & Puis fe008h . who also used a Fourier technique 
to measure rotational and macroturbulent velocities. We note, how- 
ever, that within their limited sample, the trends in both rotational 
and macroturbulent velocity mirror that of our targets. 

The only study that has sufficient overlap to allow a meaning- 
ful statistical comparison is that of Howarth et al. dl997h and in 
Fig. E] the two sets of estimates are compared for 42 stars, with 
spectral types from BO to B5. The cross-convolution method of 
Howarth et al. generally leads to larger estimates. For the slowly 
rotating stars (where our estimates are less than 70 km s _1 ), there 
is a systematic difference of approximately 30 km s _1 (with a stan- 
dard deviation of 10 km s _1 ). We ascribe this difference to the 
effect of macroturbulence, which was not explicitly included in the 
estimates of Howarth et al. For the stars rotating with intermedi- 
ate projected velocities between 70 and 140 km s , we find good 
agreement between the values obtained from the cross correlation 
and Fourier transform methods; the average difference between the 
results of the two techniques is 8 km s _1 with a standard deviation 
of 8 km s _1 . For the four most rapidly rotating stars, HD 64760, 
HD 157246, HD 99857 and HD 93827, the Fourier Tranform esti- 
mates are higher by typically 30 km s" 1 from those using the cross- 
correlation method. We would expect our estimates to be secure as 
the line profiles are now dominated by rotational broadening and 
the zeros in the Fourier Transform are well observed. It is possible 
that the difference between our values and those of Howarth et al. 
is due to gravity darkening. For rapidly rotating stars, we expect to 
see a temperature gradient acros s the star, wit h the poles being hot- 
ter than the equator (Von Zeipel ll924al . ll 924bf) . This could lead to 
Howarth et al. finding lower values for w e sin i, as their UV spectra 
are preferentially sampling the hot, slowly rotating, poles over the 
cooler, and rapidly rotating equator. 

5.2 Comparison with evolutionary models 

Our estimated projected rotational velocities are plotted against 
spectral type in Fig. [5] with the majority of stars lying on a locus 
from approximately 60 km s _1 at B0 to 30 km s _1 at B9. This de- 
crease in projected rotational velocity at later spectral type is in ac- 
cordance with evolutionary models, as discussed below. However, 
there are a number of early-B type supergiants that appear to have 
anomalously high projected rotational velocities for their spectral 
type and these will be further discussed in Sect. 15. 31 

We used a Kolmogorov-Smirnov (K-S) test to compare the 
distribution of rotational velocities in our Galactic supergiant sam- 
ple, with that found for supergiants in the SMC by Dufton et al. 
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Figure 4. Comparison of the projected rotational velocities for stars com- 
mon to both Howarth et al. 119971) and this paper. The effects of macrotubu- 
lence, which is not explicitly included in the former can be clearly seen. 

(2006). While the sample of Dufton et al. is small, with only 13 
stars, the K-S test indicates that at the 0.05 significance level, the 
distributions are the same, with a measured D statistic of 0.29. From 
this, we conclude that there is no evidence for any dependence of 
the current rotational velocities on initial metallicity, although this 
result must be qualified by again stressing the relatively small sam- 
ple sizes, and the additional uncertainty introduced by the unknown 
value of sini. 

From a theoretical perspective, the dependance of surface rota- 
tional velocity on metallic ity is not straightforward. As discussed in 
Meynet & Maeder ( 2005), both mass loss and the internal coupling 
between core and envelope affect the surface equatorial velocity. 
Lower metallicity reduces mass loss rates (and hence the loss of 
angular momentum from the star), but it also lowers the efficiency 
of the transport of angular momentum to the surface. Hence as the 
metallicity changes, there are two competing mechanisms that to 
some extent cancel each other out. 

Using the models of Meynet & Maeder d2003Ll2005l) we com- 
pared two 25 Mq models with initial velocities of 300 km s^ 1 and 
Z=0.02 and Z=0.04 (D25z20S3A and D25z40S3 respectively). The 
higher metallicity model was found to have a rotational velocity 
that was ~ 15 km s _1 lower at log T c g = 4.35, and ~ 5 km s _1 
lower at log T c g = 4. 1 . Such small differences are consistent with 
the cumulative frequency plot discussed above. The 40 Mq mod- 
els (D40z08S3, D40z20S3A and D40z40S3) show a larger spread 
of equatorial velocities for a given temperature. For example at a 
temperature of log T c ff = 4.2, the low metallicity model (Z=0.008) 
has v = 37 km s~ , the solar metallicity model has v = 43 km s , 
and the high metallicity model (Z=0.040) has v = 66 km s . Scal- 
ing these values to account for average projection, the spread in 
projected rotational velocities is ~ 20 km s , From this, we con- 
clude that while differences in metallicity may increase the scat- 
ter of measured values when compared with evolutionary models, 
these are unlikely to be large enough to be observable when com- 
paring the distribution and mea n of ro tational velocites for the SMC 
B supergiants of Dufton et al. d2006h to our Galactic sample. Fur- 
thermore, we can discount different metallicities as an explanation 
for our most rapidly rotating stars, as discussed in Section |531 

The projected rotational velocities for all our targets have 
been plotted against their effective temperature estimates (from 
Table [3j in Fig. [6] together with the predictions of the evolution- 
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ary models of Meynet & Maeder (120030 for Galactic metallici- 
ties. To aid the comparison, the observational dataset has been sub- 
divided into four initial mass ranges. The models used for the com- 
parison are D09z20S3A, D12z20S3A, D215z20S3A, D20z20S3A, 
D25z20S3A and D40z20S3A, which have Z=0.02, an initial equa- 
torial rotational velocity of 300 km s~ and masses from 9 Mo 
to 40 Mq; the theoretical rotational velocities have been scaled by 
a factor of 7r/4 to convert them to an average projected rotational 
velocity assum ing a random orientation of axes (Chandrasekhar & 
Mtinch. ll950h . 

In general there is good agreement between observation and 
theory, both in the magnitude of the projected rotational veloci- 
ties and a decrease as one moves to lower effective temperature or 
later spectral type. The agreement is particularly encouraging for 
the higher mass stars with initial mass estimates of > 20 Mq. For 
the two lower mass cohorts, the agreement is less convincing. For 
example, the projected rotational velocity estimates for our lowest 
mass sample (9-15 Mq) appear to be systematically lower than 
predicted. This discrepancy would be ameliorated if the initial stel- 
lar masses had been underestimated but as discussed in Sect. 14.51 
these estimates should be reliable to within 2 M© . Another possi- 
bility is that the projected rotational velocities have been system- 
atically underestimated. Aerts et al. d2009h simulated the effects 
of non-radial gravity-mode oscillations on the spectral line profiles 
of an early-B type supergiant and found that this could lead to an 
underestimate of the rotational velocity found by our methodology. 
Although this would again improve agreement with theory, it would 
not explain why the effect is only found in the lowest mass cohort. 
Additionally the simulations indicate that the underestimation only 
occurs in some cases and hence do not currently explain the sys- 
tematic nature of the discrepancy. It should be noted that the sim- 
ulations of Aerts et al. are for one set of atmospheric parameters 
(T e ft= 18 200 K; log g= 3.05) and an extension to other spectral 
types would be useful. Hence we conclude that supergiants with 
initial masses in the range 9-15 Mq would appear to have lower 
rotational velocities than are predicted by the evolutionary models. 
For supergiants with initial masses of 16 - 20 Mq, the agreement is 
good apart from a small number of targets that appear to be rotating 
more rapidly and which are discussed further in Sect. 15.31 

We also note that the decrease in rotational velocities seen at 
spectral type B 1 corresponds to the position of the bi-stability jump 
(Vink, de Koter & Lamers, 1999), where wind properties change 
from a fast wind with standard mass loss rates, to a slower wind 
with higher mass loss rates. It has been suggested (Vink, in prepa- 
ration) that the bistability jump may be accompanied by bi-stability 
braking, where increased mass loss at spectral type Bl removes 
sufficient angular momentum from the system to account for a dra- 
matic decrease in rotational velocities. Our results for the 20 to 25 
Mq cohort, as illustrated in Fig. [6] provide some tentative support 
for this rapid decrease in rotational velocity, with three or four tar- 
gets targets appearing to lie in this phase. 



5.3 Rapid rotators 

For the purposes of this discussion, "rapid rotators" refers to four 
targets, HD157246, HD93827, HD99857 & HD165024 that have 
larger projected rotational velocities than predicted by the evolu- 
tionary models (see Fig.(6]l. Of these, HD165024 should be con- 
sidered as a marginal case; if its mass estimate was reduced from 
17 to 12 Mofhen it would not be classed as rapidly rotating. How- 
ever, it was decided to include this in the sample of rapid rotators 
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ure 5. Projected rotational velocities as a function of spectral type. The 
additional targets that were selected as they were suspected to have a 
h projected rotational velocity have been identified. 



that would be examined in more detail for clues as to a different 
evolutionary history. 

It should be noted that these rapidly rotating objects appear to 
be relatively rare with only four identified in our complete sample. 
Indeed six targets were pre-selected based on them having anoma- 
lously large projected rotational velocities estimates for their spec- 
tral type in the analysis of Howarth et al. Excluding these stars, 
there are fifty one targets in our sample, with two having measured 
projected rotational velocities that appear to be significantly larger 
than the mean for their spectral type, viz. HD157246, HD165024. 
Thus we find that approximately 4% of our sample are anomalously 
rapid rotators. This is broadly consistent with the identification of 
only one rapid rotator in the sample of Dufton et al (2006) of 13 
SMC supergiants. Note that we may have underestimated the frac- 
tion of supergiants that are rapid rotators, as some of our sample 
could be observed pole-on. However as discussed by Hunter et al. 
d2007l) the probability of observing such objects is low due to the 
nature of the sine function. For example, only thirteen percent of 
our targets will have a value of sin i less than 0.5, assuming a ran- 
dom orientation of the rotation axis. Hence we consider that within 
the constraints of our sample size, our estimate of the fraction of 
anomalously rapid rotators among our targets is reasonable and that 
it is certainly small. 

There are at least three channels by which these stars could 
have reached their current evolutionary state: 

• they had extremely high rotational velocities on the main se- 
quence, and have spun down normally as they evolved. 

• they have not spun down as much as most other stars as they 
became supergiant stars. 

• they had a moderate velocity on the main sequence but have 
experienced an additional source of angular momentum that has 
allowed them to maintain their rotational velocity. 

The first two channels appear implausible as for the first, a 
very high main sequence rotational velocity is required, whilst for 
the second, it is unclear what mechanism would have inhibited their 
spin down. Indeed for models in this mass range the increase in 
mass loss due to the bi-stability jump naturally leads to a decrease 
in the rotational velocity. 

This leaves the possibility that these stars have an additional 
source of angular momentum. A plausible scenario would involve 
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Figure 6. Projected rotational velocities and effective temperatures for all stars, plotted against the evolutionary tracks of Meynet and Mader 120031) . Note that 
the rotational velocities from the evolutionary tracks have been scaled by n/A to account for projection, apart from the red dotted line in the second panel that 
displays the evolution of a star with an initial mass of 15 Mq . Stars (which have had their masses estimated as described in Section 14. 5\ have been binned 
accordingly into four mass ranges 9M Q < M < 15M , 16M < M < 2OM , 21M < M < 25M and 26M < M < 4OM . These mass bins were 
chosen so as to facilitate comparison with the available evolutionary models. 



binarity and mass transfe r (Meibon, Mathieu & Srassun, 120071 ; 
Langer et al. I2008L 120031 ; Koenigsberger et al. l2003h . The SIM- 
BAD online astronomical database was checked to see if any of 
the four rapid rotators were listed as spectroscopic binaries (see 
Fig|3j. None of the targets were so classified, and in the absence 
of the time-resolved spectra that could indicate binarity from ra- 
dial velocity variations, the spectral lines in the rapid rotators were 
re-examined but no evidence was found for any asymmetry that 
might imply the presence of a binary companion. This is not in it- 
self evidence that these are single stars, as the luminous nature of 
B-supergiants would make the detection of fainter binary compan- 
ions difficult. 



In addition to these three scenarios, we must also consider the 
possibility that the models used are simply not appropriate for these 
objects. It would be surprising if this was the case however, as the 
agreement between theory and observation is excellent for most of 
our sample, and there is nothing particularly unusual about the ob- 
served properties of our rapid rotators (besides their large v a sin i) 
that would lead us to suspect the models are not applicable. 



5.4 Velocity fields 

The atmospheres of early-type supergiants appear to contain sig- 
nificant velocity fields. These are often characterised as microtur- 
bulence (estimated in Sect. [43} and macroturbulence (estimated in 
Sect. 13 . 21 ) . The former is incorporated into the line opacity profile 
and hence directly influences the strength of the metal line, whilst 
the latter is normally modelled by including additional broadening 
via a convolution of the theoretical line profile. Assuming that these 
effects are indeed caused by velocity fields, they represent the ex- 
tremes of velocity fields with different scale lengths in the stellar 
atmosphere. 

The microturbulent velocities deduced for our sample range 
from approximately 12 to 24 km s _1 (with two cases where no reli- 
able estimate was possible). There is some evidence that the micro- 
turbulent velocity may be anti-correlated with the effective temper- 
ature but this result should be treated with caution for at least two 
reasons. Even with very high quality data, there are significant ran- 
dom errors associated with the estimation of the microturbulence. 
For example, the estimate of the microturbulence for HD 150898 is 
20 km s _1 . However arbitrarily decreasing the equivalent width of 
Sim line at 4552 A by 5 per cent, whilst increasing that for the line 
at 4575 A by 5 per cent would decrease the estimate to 15 km s _1 . 
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Additionally the microturbulence estimate can vary depending on 
the ionic species an d the set of lines considered (s ee, fo r exam- 
ple, Vranck en et al. 1 1991 l200d: Simon-Diaz et al. 120061 : Hunter 
et al. 120071) . Indeed Vrancken et al. reported differences of up to 
9 km s _1 in the estimates obtained from different species. Addi- 
tionally the strength of the Si III triplet varies with spectral type 
leading to different degrees of saturation and hence of sensitivity to 
the microturbulence parameter. 

Recently Cantiello et al. (2009) have discussed the obeser- 
vational consequences of sub-surface convection zones in massive 
stars. These zones arise from peaks in the opacity (as a function of 
temperature) that arise from the ionization of helium and iron. The 
latter is found to be more prominent at higher luminosity and lower 
gravities corresponding in the HR-diagram to the region occupied 
by O and B-type superg iants. U sing the analyses of Trundle et al. 
d2007t) and Hunter et al. d2008ij) of B-type stars observed in the ESO 
VLT-FLAMES survey of massive stars (Evans et al. l2005h . Cantiello 
et al. find a correlation between the magnitudes of the convection 
velocities implied by their models and the estimated microturbu- 
lent velocities. In particular objects with microturbulent velocities 
of more than 10 km s _1 are in a region of the HR-diagram where 
the mean convective velocity within 1.5 pressure scale heights of 
the upper boundary of the iron convection zone (< v c >) is greater 
than 2.5 km s . 

All our targets lie in the region of the HR-digram where the 
convection due to the iron opacity is pronounced. Indeed as well 
as lying within the contour for < v c >= 2.5 km s -1 , they all lie 
within the contour for 3.75 km s _1 which was the largest mean 
convection velocity regime considered by Cantiello et al. For their 
Galactic sample, Cantiello et al. only considered nine stars with five 
of them having estimates of the microturbulent velocity in excess 
of 10 km s . For our sample, we have been able to estimate the 
microturbulent velocity in 54 (out of 57) targets. The large values 
deduced coupled with their position within the HR-diagram sup- 
ports the conclusion of Cantiello et al. that there is a physical con- 
nection between sub-photospheric convection and the small scale 
stochastic velocity fields in early-type stars. 

As discussed in Sect. 13.21 a macroturbulent velocity was in- 
troduced in order to allow for the excess broadening in the ob- 
servational metal line profiles as has been found in other studies 
(see, for example, Howarth et al. 1 19971 ; Ryans et al. 2003). The 
velocity field was assumed to have a Gaussian distribution and in 
general this appeared to be consistent with the observations. Es- 
timates were obtained for all but six of our targets. In four cases 
(HD 64760, HD 152667, HD 157246, HD 93827) the line profile is 
dominated by rotational broadening making the estimation of the 
relatively small macroturbulent component unreliable. One target, 
HD 83183, appeared to have a very small degree of macroturbulent 
broadening, possibly due to its having a luminosity class II. For the 
sixth target, HD 77581, the fits to the spectra were relatively poor 
and a reliable estimate was therefore not possible. 

The errors for the macroturbulence quoted in Table [3] are the 
standard deviations of the estimates from individual features. As 
such they do not include systematic errors and implicitly assume 
that there is a unique value for the macroturbulence. A potentially 
significant source of systematic error is the choice of the intrin- 
sic absorption line profiles. These were taken from spectra calcu- 
lated from our grid of non-LTE TLUSTY model atmospheres as de- 
scribed in Sect. 13.21 Tests using theoretical spectra from different 
grid points indicated that uncertainties in the values adopted for 
the effective temperature or surface gravity were unlikely to be sig- 
nificant. However the choice of the microturbulence (the TLUSTY 
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peratures estimates for all stars in our sample for which reliable estimates 
were o btained (marked +). Also shown are the estimates of Dufton et al. 
1 2006) for Galactic supergiants (marked x) and SMC supergiants (marked 
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grid was calculated for values of 0, 5, 10, 20 and 30 km s _1 ) was 
more important as it is this parameter that effectively determines 
the shape of the theoretical line profile. Again tests showed that 
for the early-B type stars, the value adopted for the microturbu- 
lence had little effect of the macroturbulence estimates. This was 
because the former had values (10-20 km s _1 ) that were signifi- 
cantly smaller than the latter (typically 60 km s" 1 ). In such cases 
the intrinsic profile is relatively narrow compared with the macro- 
turbulent broadening and hence its choice is not critical. However 
for late B spectral types, the magnitudes of the microturbulence and 
macroturbulence become comparable and in this case the choice of 
the intrinsic line profile becomes more important. 

For example, the grid point adopted for HD51309 had a mi- 
croturbulence of 20 km s _1 being the nearest value to the obser- 
vational estimate of 18 km s~ . If instead we had adopted a mi- 
croturbulence of 10 km s - , our estimate for the macroturbulence 
would have increased from 20±3 km s _1 to 31±5 km s . This 
result is not surprising as we are modelling the total broadening of 
the lines using micro- and macroturbulent velocity fields (together 
with a rotational component). A decrease or increase in magnitude 
of one field will inevitable anti-correlate with the estimate obtained 
for the other field. However the uncertainty in the macroturbulence 
estimates for the later B spectral types may be larger than that im- 
plied by the statistical error estimates quoted in Table [3] 

The estimates of the macroturbulent velocity show a decrease 
from approximately 70 km s _1 at a spectral type of B0 to approxi- 
mately 20 km s _1 at a spectral type of B8. Two targets (HD 52089 
and HD 159100) appear to have anomalously low estimates. How- 
ever these estimates will depend on the adopted microturbulences 
as discussed above. Additionally both these objects have relatively 
large surface gravities for their effective temperature, together with 
lower estimated masses than most of our sample; hence their lower 
estimates may reflect their different evolutionary status. The de- 
crease of macroturbulence with effec tive temperature has been 
found previously by Ryans et al. ( 2002) for 1 1 Galactic supergiants 
and by Dufton et al. (2006) for 13 SMC supergiants. These results 
have also been plotted on the Fig. [7] and in general the estimates 
obtained from all three st udies a re consistent. 

Recently Aerts et al. d2009h have provided a physical explana- 
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tion for this excess broadening in terms of the non-radial gravity- 
mode oscillations in the stellar atmosphere. They simulate the ef- 
fect of these oscillations on the profile of the Si III line at 4552 
A in a B2 la type supergiant (T cff ■= 18 200 K; log g= 3.05). Addi- 
tional broadening is found which is modelled by rotational broad- 
ening coupled with a Gaussian macroturbulent velocity field. For 
the former, the estimates can be lower than those of the models 
(as discussed in Sect. 15. 2\ , whilst estimates of the macroturbulence 
velocity ranging from km s _1 to approximately 50 km s _1 are 
found. This range is consistent with our observational estimates in 
Fig. Hand additionally the range of values found in the simulations 
(due to the time dependent contribution of the different modes) may 
also at least partially explain the scatter found in our estimates, in- 
cluding the anomalously low estimates found for HD 52089 and 
HD 159100. 

5.5 Nitrogen abundances 

Rotation is generally considered to be crucial for developing theo- 
retical models of massive star evolution as it can induce mixing of 
nucleosynthetically processed material from the stellar core into the 
photosp here ( see, for example, Heger & Langer 2000; Meynet & 
Maeder l200d) . Such models have been used to explain the ratio of 
blue to red supergiants (Maeder & Meynet 2001) and Wolf-Rayet 
populations at different metallicities (Meynet & Maeder 2005]; Vink 
& de Koter l2005h . A natural consequence of such mixing is an en- 
hancement of nitrogen (together with smaller changes in the car- 
bon, oxygen and helium abundances) at the stellar surface. Hence 
nitrogen abundances are often used as an indicator of the amount of 
mixing, and in particular rotationall y-indu ced mixing, that has oc- 
curred (see, for example, Brott et al. l2008l) . Howeve r recent results 
from the VLT-FLAMES survey (Evans et al. l2005h have indicated 
that other mechanisms may affect m ixing of mate rial between the 
interior and the surface (Hunter et al. l2008bl . 120091) . 

We have estimated the surface nitrogen enhancement for our 
targets by comparing our nitrgen abundances with those found for 
Galactic B-ty pe star s in the clusters NGC 3293 and NGC 4755 
(Trundle et al. l2007h . This study was choosen as it used the same 
grid of model atmospheres and methodology to that adopted here. 
To ensure consistency with our own results, we have only used their 
abundance estimates from the 3995 A line, whilst to minimise the 
effects of mixing, we have only included stars of luminosity class 
V-III. A mean abundance of 7.47±0.14 dex was obtained from a 
sample of 15 stars. This is slightly smaller than the mean abundance 
of 7.59±0. 19 dex found when all the nitrogen lines are considered 
but is consistent with the comparison of nitrogen abundances given 
in Table [5] In Fig. [8] the surface nitrogen abundance for our tar- 
gets is plotted against projected rotational velocity, together with 
the base abundance deduced from the results of Trundle et al. 

There would appear to be no correlation of nitrogen abun- 
dance with current projected rotational velocity. This is not surpris- 
ing given the ambiguity arising from the unknown angle of incli- 
nation and more importantly the large change in rotational velocity 
that occurs when a star leaves the main sequence (as is apparant, 
for example, in the theoretical models shown in Fig. [6). We have 
attempted to allow for the latter effect by estimating the projected 
rotational velocity of the progenitor main sequence object. We have 
assumed that all the main se quence precursors have a typical sur- 
face gravity of 4.2 dex (Cox (2000h ; we note that the actual value 
adopted is not critical to our methodology. We have then assumed 
that the decrease in rotational velocity as a star leaves the main 
sequence is due solely to its increase in radius. Hence we are as- 



suming that the effects of mass loss in both decreasing the stellar 
moment of inertia and in removing angular momentum are neg- 
ligible. In such circumstances the radius will simply scale as the 
square root of the surface gravity and we can estimate the ratio of 
the main sequence to the current radius. Additionally, if we assume 
solid body rotation and that the stellar moment of inertia scales as 
the radius squared, we can estimate the ratio of the main sequence 
to the current projected rotational velocities and hence deduce esti- 
mates for the former. 

To test the validity of our approach, we considered the de- 
crease in equatorial rotational velocity as the radius of a star in- 
creases, as predicted by the Geneva stellar evolutionary models. 
These tests indicated that our simplification was reasonable, with 
best agreement for the 20Mq model. For this model, evolving from 
the main sequence (defined here as when the star has a surface grav- 
ity of 4. 1 dex) to the supergiant phase when the star has a surface 
gravity of 2.9 dex, its rotational velocity decreases by a factor of 
3.9, whilst its radius increases by a factor of 3.8. For more massive 
models, the agreement is less satisfactory, with a 25 Mq model 
increasing in radius by a factor of 3.8 as it moves from the main 
sequence to the supergiant stage, but decreasing its rotational ve- 
locity by a factor of 5.6. At lower masses the agreement is poorer, 
with a 15 Mq model increasing its radius by a factor of 4, while 
decreasing its rotational velocity by a factor of 2.7. These discrep- 
ancies are most likely due to mass loss, which we have ignored in 
our simplification. Furthermore, we have assumed solid body ro- 
tation, whereas in the more sophisticated Geneva models, the cou- 
pling between core and envelope and the changing size of the core 
will affect the transport of angular momentum within the star, and 
hence rotational velocity. 

In Fig. [9] our nitrogen abundances are plotted against these 
main sequence projected rotational velocities (hencefore desig- 
nated as "scaled projected rotational velocities"). We have excluded 
the four targets that we believe may not have been undergone nor- 
mal single star evolution (see Sect. 15. 3t but our conclusions would 
not be significantly affected if they were included. We also ex- 
cluded HD 152236 as its surface gravity could not be measured, 
and HD 64760 as its scaled rotational velocity is unrealistically 
large. We emphasize that the scaled projected rotational velocities 
must be considered as subject to significant uncertainties given the 
approximations made. Additionally they do not remove the ambi- 
guity introduced by the unknown angle of inclination. However we 
believe that they are useful as the large increase in radius is likely 
to be an important mechanism as the star evolves from the main 
sequence. 

A correlation is now seen between the scaled projected rota- 
tional velocity and nitrogen abundance. A linear least squares fit 
(shown in Fig. [6] as a solid line) implies an increase in nitrogen 
abundances of approximately 0.6 dex (i.e. a factor of four) over our 
range of scaled projected rotational velocities. The slope of the fit is 
found to be (1.1 ±0.3) x 10 -3 dex / km s _1 , which is not consistent 
with a null result. Additionally the asymptotic nitrogen abundance 
at low velocities of approximately 7.7 dex is in reasonable agree- 
ment with our baseline nitrogen abundance of 7.47 dex. Indeed as 
some of the stars with low projected rotational velocities will have 
low angles of inclination (and hence relatively large projected ve- 
locities) a small positive offset from the baseline abundance would 
be expected. We used a Monte-Carlo technique to try and test the 
validity of the correlation. The value of v c sin i was randomly var- 
ied within one standard deviation of the mean for each star, and 
nitrogen abundance and surface gravity varied by 0.2 and 0.1 dex 
respectively, and a least squares fit made to the data. For all iter- 
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Figure 8. Surface nitrogen abundaces versus projected rotational velocity. 
Solid line is base abundance for early B stars while on the main sequence, 
based on VLT-FLAMES data as discussed in text. 
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Figure 9. Surface nitrogen abundaces versus estimated main sequence pro- 
jected rotational velocity. Also shown is a linear least squares fit. 



ations, a positive correlation between scaled v e sin i and nitrogen 
abundance was observed, with a nitrogen enhancement of between 
0.3 and 0.6 dex. Based on these results, and within the limits of our 
approximations, there would appear to be a statistically significant 
correlation between inferred main sequence rotational velocity and 
surface nitrogen abundance at the supergiant stage. 

These results are consistent wit h stell ar evolutionary calcula- 
tions. For example, Heger & Langer (2000) discussed the evolution 
of a star with an initial mass of 20 Mq and an equatorial velocity 
of 206 km s~ 1 . They find an increase of a factor of two or more in 
the surface nitrogen abundance at the end of core hydrogen burn- 
ing. T he Geneva group stellar evolution models (Meynet & Maeder 
120001) also find similar nitrogen enhancements at the surface of 
rapidly rotating massive stars. For example, a 20 Mq model ro- 
tating with an initial velocity of 200 km s _1 , leads to an increase in 
the nitrogen to carbon ratio at the surface by a factor of four; when 
the initial rotational velocity is increased to 300 km s _1 , this in- 
creases to a factor of approximately seven with most of the change 
arising from an enhanced nitrogen abundance. Given the uncertain- 
ties in the initial rotational velocities of our objects, we conclude 
that there is satisfactory agreement between our results and models 
of rotational mixing. However our results cannot exclude the types 
of inconsistencies found by Hunter et al. d2008bh in their study of 
LMC B-type stars. 

Although we have compared our results with single star evolu- 
tionary models, we note that in some cases the surface composition 
may also have been affected by binary star evolution. For exam- 
ple, while attempting to model the binary progenitor of SN 1993J, 
Stancliffe & Elridge d2009h found that a blue supergiant in a bi- 
nary system could undergo a surface nitrogen enrichment of up to 
a factor of ten. 

5.6 Differential Rotation 

The relatively simple Fourier transform analysis employed in this 
study is based on the assumption that stars rotate as a solid body. 
While observations of o ur ow n sun have long demonstrated this to 
be false (eg. Schou et al. ll998h . it nonetheless remains an appropri- 
ate assumption when confronted with sp ectra o f limited resolution 
and signal-to-noise. Reiners & Schmitt d2002h have modeled the 
dependence of the first (qi) and second (q 2 ) zeroes in the Fourier 



transform of rotationally broadened spectral lines, and have estab- 
lished criteria for solar and anti-solar like differential rotation. Us- 
ing this technique, differential rotation has been identified in A-type 
main sequence stars by Reiners & Royer ( I2003h . 

Differential rotation has not, however, previously been ob- 
served in stars of comparable spectral type and luminosity class to 
our sample. We have therefore attempted to identify the signatures 
of differential rotation in our sample. Only the fast rotating stars 
were considered, as a large projected rotational velocity allows the 
reliable identification of more than one zero in the Fourier Trans- 
form. Furthermore, the angle of inclination of the most rapidly ro- 
tating stars can be constrained, as v ~ v c sin i when i ~ 90° . This 
is important as pole-on rotation may mimic the characteristics of 
differential rotation in the Fourier transform of absorption lines. 

For the majority of the targets, the second zero was indistin- 
guishable from the high frequency noise in the transform. However, 
for five targets (HD152667 (q 2 /qi = 1.77), HD64760 (q 2 /qi = 
1.75), HD139518 (q 2 /qi = 1.81), HD100841 (q 2 /qi = 1.74) 
and HD157246 (q 2 /qi = 1.87 ± 0.08)), the position of the 
second zero was measurable. Reiners & Schmitt give a value of 
92/51 > 1-83 as being indicative of anti solar-like differential ro- 
tation (where the poles rotate faster than the equator), and a value 
of q 2 /qi < 1-72 for solar-like differential rotation. Hence of these 
stars, only HD 157246 satisfied the criteria for differential rotation, 
and in Table [6] we summarize the measurements for this star. It is 
important to note however, that the value of q 2 / q\ determined for 
HD 157246 is within a standard deviation of the threshold for dif- 
ferential rotation. 

While anti solar-like differential ro tation has recently been ob- 
served by, for example, Strassmeir et al. (2003), from Doppler maps 
of a K2 giant, it is unclear what the physical mechanism behind 
this is. It has been proposed that meridional circulation, driven by 
a strong temperature gradient, could lead to anti solar-like differen- 
tial rotation (Kichatinov & Riidiger, 2005). 



6 CONCLUSION 

From a sample of more than fifty early type supergiant stars, we 
find a range of rotational velocities that in general are in good 
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Table 6. Details of lines used for analysis of differential rotation in 
HD157246 



Species 


A (A) 


v e sin i 


?l 


92 


32. 
91 


He I 


3819.61 


264 


0.194 


0.359 


1.85 


Nil 


3995.00 


279 


0.177 


0.335 


1.89 


He I 


4026.19 


269 


0.181 


0.346 


1.91 


CII 


4267.00 / .26 


255 


0.181 


0.335 


1.85 


Mgll 


4481. 13/. 33 


255 


0.171 


0.296 


1.73 


Sim 


4552.62 


277 


0.156 


0.296 


1.90 


on 


4661.63 


279* 


0.150 


0.301 


2.01 


He I 


4921.93 


271 


0.148 


0.274 


1.85 



* Poor quality line, hence q\ and q2 are uncertain. 



agreement with evolutionary models. For four stars, we find a sig- 
nificant difference between their measured values and the expected 
rotational velocities for their current evolutionary status. These dis- 
crepancies cannot be explained in terms of projection effects and 
are unlikely to be caused by observational uncertainties, and we 
conclude that the most likely explanation for their behavior is bi- 
narity. In addition to the projected rotational velocities, we have 
measured the macroturbulence and microturbulence in most of our 
targets. For the former, we find good agreement with the predic- 
tions of Cantiello et al. for convection driven by the opacity 
of iron group elements. The latter may be a manifestation of the 
large number of non-radial gravity-mode stellar oscillations (Aerts 
et al. 2009) present in the stellar photosphere and our results are 
consistent with their prediction for a B2 la supergiant. 

We find a range of nitrogen abundances stetching from the 
Galactic B-type baseline abundance to an enhancement of more 
than 1 dex. No correlation of abundance is found with the current 
projected rotational velocity. However a correlation is found with 
the estimated main sequence projected rotational velocity and this 
is consistent with current evolutionary models for single rotating 
stars. 

We have attempted to determine if there was any evidence of 
differential rotation in our sample, and find marginal evidence that 
one of our most rapidly rotating targets is undergoing anti-solar like 
differential rotation. 
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